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We prove exponential decay of the off-diagonal correlation function in the two-dimensional homo- 
geneous Bose gas when a^p is small and the temperature T satisfies 

T> . 
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Here, a is the scattering length of the repulsive interaction potential and p is the density. To leading 
order in a^p, this bound agrees with the expected critical temperature for superfluidity. In the 
three-dimensional Bose gas, exponential decay is proved when 
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> 5\/ap^ , 



where ri"' is the critical temperature of the ideal gas. While this condition is not expected to be 
sharp, it gives a rigorous upper bound on the critical temperature for Bose-Einstein condensation. 
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I. INTRODUCTION 

Quantum many-body effects due to particle interac- 
tions and quantum statistics make the Bose gas a fasci- 
nating system and a challenge to theoretical physics. It is 
increasingly relevant to experimental physics, especially 
after the first realization of Bose-Einstein condensation 
in cold atomic gases It displays a stunning physical 
phenomenon: superfluidity. Several mechanisms that are 
present in the Bose gas also play a role in interacting elec- 
tronic systems and in quantum optics. 

Both the two-dimensional and the three-dimensional 
gas have physical relevance, and they behave rather dif- 
ferently. We consider them separately here. Throughout 
the paper, we shall assume that units are chosen in such 
a way that h ~ 2m = fee = b where m is the particle 
mass. 



A. The two-dimensional Bose gas 

There is no Bose-Einstein condensation in the two- 
dimensional Bose gas at positive temperature, as was 
proved by Hohenberg more than forty years ago.^Tn con- 
trast to higher dimensions, the ideal Bose gas offers no 
intriguing features in two dimensions. But the interact- 
ing gas is expected to display a Kosterlitz-Thouless type 
transition from a normal fluid to a superfluid, where the 
decay of off-diagonal correlations goes from exponential 
to power law. The critical temperature Tc depends on 
the scattering length a of the interaction potential, which 
we consider to be repulsive. For dilute gases, i.e. when 



a'^p <C 1, Popo\0 performed diagrammatic expansions in 
a functional integral approach, flnding that 



ln|ln(a^p)| 

This formula was conflrmed by Fisher and Hohenberg 
using Bogoliubov's theory, and by Pilati et elW using 
Monte-Carlo simulations. No rigorous proof is available 
to this date, however. 

In this article we prove in a mathematically rigorous 
fashion that there is exponential decay of the off-diagonal 
correlation function when the temperature satisfles 

^ ln|ln(a2p)| ^ *^ ( ln|ln(a2p)| )) 

for small a^p. Thus we prove that Tc cannot be big- 
ger than the conjectured value Q, to leading order in 
a^p. The main novel ingredient in our proof is a rigorous 
bound on the grand-canonical density of the interacting 
Bose gas. This is explained in the next section. 

B. The three-dimensional Bose gas 

A three-dimensional Bose gas is interesting even in the 
absence of particle interactions. Bose-Einstein conden- 
sation takes place at the critical temperature T^^'' — 
47r(p/C(|))^/^ (where C(f ) « 2.612, with C the Riemann 
zeta function) . The effects of particle interactions on the 
critical temperature have been studied by many authors. 
A consensus has been reached in recent years but it is 
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tenuous; we give a survey of the main results, both for 
historical perspective and in order to gain a sense of the 
solidity of the consensus. Let ATc = Tc - Ti°^ denote 
the change of the critical temperature. 

1953 FeynmarP argued that interactions increase the ef- 
fective mass of the particles and hence decrease Tc , 
i.e, ATc < 0. 

1958 Lee and Yan^' predict that the change of criti- 
cal temperature is linear in the scattering length, 
namely 



ATJT^ 



cap 



1/3 



No information on the constant c is provided, not 
even its sign. 

1960 Glassgold, Kaufman, and WatsoiP find that the 

critical temperature increases as 
C(api/3)i/2 with C> 0. 



1964 HuangiSl gives an argument suggesting that 
ATe/Ti°^ « C(api/3)3/2 ^i^l^ 

1971 A Hartree-Fock computation shows that ATc < 
(Fetter and WalecksPl). 

1982 A loop expansion of the quantum field representa- 
tion gives ATe/Ti"^ w -3.5(api/3)i/2 (ToyodcP). 

1992 By studying the evolution of the interacting Bose 
gas, Stoof^"^ finds that the change of critical tem- 
perature is linear in the scattering length with 
c= 167r/3C(3/2)4/3 = 4.66. 

1996 A diagrammatic expansion in the renormalization 
group yields ATc > (Bijlsma and Stoof^**). 

1997 A path integral Monte-Carlo simulation yields c = 
0.34 ± 0.06 (Griiter, Ceperley, and Lalo^. 

1999 A virial expansion leads to c = 0.7 (Holzmann, 
Griiter, and Lalo^^l). Another virial expan- 
sion leads Huan#3 to conclude that ATc/Ti"^ w 
3.5(ap^/'^)^/^. Interchanging the limit a ^ 
with the thermodynamic limit, and using Monte- 
Carlo simulations, Holzmann and KrautlP^ find 
c = 2.3 ± 0.25. The dilute Bose gas can be mapped 
onto a classical field lattice model (Baym et. al.'^; 
a self-consistent approach then yields c = 2.9. 

2000 An experimental realization by Reppy et. al.l^ 
yields c = 5.1 ± 0.9. It was later pointed out that 
the estimation of the scattering length between par- 
ticles was not correct, however. ^-'^ 

2001 Arnold and Moore,'^^' and Kashurnikov, Prokof'ev, 
and Svistuno\^^ performed numerical simulations 
on the equivalent classical field modelfl^l the former 
get c = 1.32±0.02 and the latter get c = 1.29±0.05. 



2003 A variational perturbation theory performed by 
Kleinert'23' yields c = 1.14 ± 0.11. 

2004 By studying the classical field model with varia- 
tional perturbations, KasteninJ^ll gnds c = 1.27 ± 
0.11. A path integral Monte-Carlo simulation by 
Nho and LandaiPa yields c = 1.32 ± 0.14. 

The last articles essentially agree with one another, and 
also with more recent articles .El The case for a linear cor- 
rection with constant c « 1.3 is made rather convinc- 
ingly; it is not beyond reasonable doubt, though. Notice 
that the constant c is universal in the sense that it does 
not depend on such special features as the mass of the 
particles or the details of the interactions. (The mass 
enters the scattering length a, however.) 

The question of the critical temperature for interact- 
ing Bose gases is reviewed in Baym et. alW^ and in 
Blaizot.*^ A comprehensive survey on many aspects of 
bosonic systems has been written by Bloch, Dalibard, 
and Zwerger.l^fil This question is also mentioned in addi- 
tional articles dealing with certain perturbation methods. 
The value of c is assumed to be known and its calculation 
serves to test the method. Some of these references can 
be found in Blaizot.^^^ 

In this article we give a partial rigorous justification of 
the results in the literature by proving that off-diagonal 
correlations decay exponentially when 



T-T 



(0) 



T 



(0) 



> 5.09v/api/3 (l + 0(v/ap^)) • (3) 



In particular, there is no Bose-Einstcin condensation 
when ([3]) is satisfied. This rigorous result is not sharp 
enough to disprove any of the previous claims that have 
been just reviewed, although it gets close to Huang's 1999 
result. As in the two-dimensional case, the proof is based 
on bounds of the grand-canonical density for the inter- 
acting gas. 



C. Outline of this article 

In the next section, we shall explain how the exponen- 
tial decay of correlations can be deduced from appropri- 
ate lower bounds on the particle density in the grand- 
canonical ensemble. These bounds will be proved in the 
remaining sections. In Section |III[ we shall state our 



main result. Theorem III.l and we shall explain the pre- 
cise assumptions on the interparticle interactions under 
which it holds. Our main tool is a path integral represen- 
tation which is explained in detail in Section [TV] Finally, 
in Section |V] we investigate certain integrals of the differ- 
ence between the heat kernel of the Laplacian with and 
without potential, and obtain bounds that are needed to 



complete the proof of Theorem III.l 
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II. DECAY OF CORRELATIONS 



We consider the grand-canonical ensemble at chemi- 
cal potential ^ and we denote the fugacity hy z — e^^^ . 
Let 7(0;, y) = {a^ {x)a{y)) denote the reduced one-particle 
density matrix of the interacting system, and 7^"-' the one 
of the ideal gas. An important fact is that, when the in- 
teractions are repulsive, we have 



7(0;,?;) ^ -i''°\x,y) 



(4) 



for any < 2; < 1. See Bratteli-Robinson,'^ Theo- 
rem 6.3.17. In d spatial dimensions. 



v(0) 



{x,y) = ^ 



which behaves like exp(— -y/— In z \x — y\) for large 
|a; — ?/|. That is, off-diagonal correlations decay exponen- 
tially fast when z < I. In particular, the critical fugacity 
satisfies Zc ^ 1. 

Next, let p{z) denote the grand-canonical density of 
the interacting system (it depends on /? as well, although 
the notation does not show it explicitly), and let 



p(°)(z) = (47r/3)-'^/2g,/2(z) 



(5) 



the density of the ideal system. Here, the function 3^/2 
is defined by 



9~ 



n > 1 



(6) 



The density p{z) is increasing in z. Then a sufficient 
condition for the exponential decay of correlations is that, 
for some z < 1, 



P < P{z) 



(7) 



The obvious problem with this condition is that the den- 
sity p{z) for the interacting system is not given by an 
explicit function. Our way out is to obtain bounds for 
p{z) (see Theorem III.l below) and to use them with 
z < 1 suitably chosen. 



A. Two dimensions 

We now explain the proof of exponential decay of corre- 
lations under the condition (|2]) for = 2. We show below 
(see Theorem III.l and the following remarks) that the 
density satisfies the lower bound 



p{z) ^ p(«)(z)- 



C |ln(l-z)| 



1 



47r/3 



|ln(a2//?)| 



(8) 



for some constant C > and for a(3^^^^ small enough. 
Here, a denotes the two-dimensional scattering length, 
which can be defined similarly to the three-dimensional 



case via the solution of the zero-energy scattering 
equation .1221311 

In two dimensions, p(°''(z) — —(AttP)^^ ln{l — z). For 
the choice z = zq with 



zo = 1 



ln|ln(aV/3)| 



ln(a7/3)| ' 
the criterion ^ is fulfilled when 

^ ^ ln|ln(aV/3)| ^^ ^^Inln | ln(aV/3)| 



47r/3 



ln|ln(aV/3)| 



Since [3 ~ 1/T, one can check that this is equivalent to 
the condition 

The situation is illustrated in Fig. [T| with qualitative 
graphs of p^^\z) and p(z). The critical fugacity z^ is 
known to be larger than 1. Our density bound holds for 
z < 1, and this yields the lower bound po for the critical 
density. It turns out to be equal to the conjectured crit- 
ical density (determined by Eq. ([T])) to leading order in 
the small parameter a/3~^/^. 




FIG. 1: Qualitative graphs of the grand-canonical density for 
d = 2. The shaded area represents our lower bound for the 
interacting density — the darker area is the function defined 
in Eq. Q and it extends to the lighter area by monotonicity 
of the grand-canonical density. Our lower bound po for the 
critical density is obtained by choosing zq = 1 — ■ 



B. Three dimensions 

We shall prove exponential decay under the condition 
([3]), where the constant 5.09 is really 



27/2^1/2 

3C(3/2)7/6 



^23/2 + C(3/2) 



5.09. 
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It is more convenient to consider the change in the critical 
density rather than in the temperature. Inequahty (|3| is 
equivalent to 



(0) 

.(0) 



-A'v/a^FV2 (l + 0(v/a/FV2)) , (9) 



where pi^'^ — is the critical density of the ideal 

Bose gas at temperature T, and where the constants A 
and A' are related by 

3 C(3/2)V'^ 
^ -2 (47r)i/4 ^-4.75. 

We show below that the lower bound 



{27: py 



(2^/2 + c(i)) 



— In ; 



C 



(10) 



holds for some positive constant C and a/3 small 



enough (see Theorem III.l and the following remarks). 

i 

(11) 



We use d(73/2/dz = z ^1/2(2), as well as the bound 
51/2(2:) <C 

to obtain 



Vi 



dt 



- In z 



JO) 



TT ds 



— In s s 
{ATr)-^p-'^\/~lnz. 



The criterion ^ is thus fulfilled when 
p < p(°)(l)-(4^)-i/3-5\/^lir^ 



-In; 



C 



for some z < 1. The right side of this expression depends 
on z only through w = \/~ In -z. Since the minimum 
of Aw + ^ over w > is 2VAB, we get the condition 
([9]). Notice that the optimal choice of z is zq = 1 ^ 
^-1/2(23/2 _^ (^(3/2)) a/3- to leading order in a/3-1/2. 

The three-dimensional situation is illustrated in Fig. [2] 
The critical fugacity Zc is larger than 1 but our density 
bound holds for z < 1. Our lower bound for the criti- 
cal density, po, is close to the conjectured expression for 
small 0/3--^/^^ 



III. RIGOROUS DENSITY BOUNDS 

We are left with proving the lower bounds ([8| and 
These will be an immediate conse- 



(10 1, respectively, 
quence of Theorem |III. 1| below. In order to state our re- 
sults precisely, we shall first give a definition of the model 
and specify the assumptions on the interaction potential. 



l/2„-7/4 

~ a p 




FIG. 2: Qualitative graphs of the grand-canonical density for 
d = 3. The shaded area represents our lower bound for the 
interac ting density — the darker area is the function defined 
in Eq. ( 10 1 and it extends to the lighter area by monotonicity 
of the grand-canonical density. Our lower bound po for the 
critical density is obtained by choosing zo = I — Caf3~^^^. 



The difference between p, 
order af3~^. 



(0) 



and pc is expected to be of the 



This makes it necessary to adopt a precise mathemati- 
cal tone from now on. We do so in order to make the 
results accessible also to readers with a more mathemat- 
ical background. 

Let A C M'* be an open and bounded domain. The 
state space for N bosons in A is the Hilbert space 
Lgyj^(A^) of square-integrable complex- valued functions 
that are symmetric with respect to their arguments. The 
Hamiltonian is 



N 



Ha. 



N 



U{xi-Xj), 

1 ^ i<j <: N 



with Ai the Laplacian for the i-th variable, with Dirichlet 
boundary conditions on the boundary of A. The repul- 
sive interaction is given by the multiplication operator 
U{x) ^ 0. We assume that U is radial and has finite 
range, i.e., [/(a;) = for |a;| > i?o- No regularity is as- 
sumed, however; we only require that the Hamiltonian 
defines a self-adjoint operator on an appropriate domain, 
and that the Feynman-Kac formula for the heat kernel 
applies. In particular, U is allowed to have a hard core. 
The scattering length of U is denoted by a. 
The grand-canonical partition function is 



Z = Z(/3,A,z)- J2 



z^Tr e-'^^*-" 



N ^0 

The thermodynamic pressure is defined by 

1 



p(/3,z) 



/3|A| 



lnZ(/3,A,z), 
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and the density is given by 



p{z) = (3z—p{(3,z). 



(12) 



We always work in finite volume A. The existence of 
the thermodynamic limit for the pressure, density and 
reduced density matrix is far from trivial. In particular, 
the limit for the latter has only been proved when z is 
small enough. ■^^ This is of no relevance to the present arti- 
cle, however, since our bounds apply to all finite domains 
uniformly in the volume. The one-particle reduced den- 
sity matrix can be written in terms of the integral kernels 
of the operators e~^^'^-" as 

7(2;, y) = ^ Nz^ I dx2--- dxjv 



N > 1 



X e 



{X,X2, ■ . .,XN;y,X2, ■ ■ .,Xn) 



Relatively few rigorous results on interacting homo- 
geneous Bose gases are available to this date. The 
only proof of occurrence of Bose-Einstein condensation 
deals with the hard-core lattice model at half-filling.'22inSl 
Roepstorfl^ used Bogohubov's inequality to get an up- 
per bound on the condensate density. Several as- 
pects of Bogoliubov's theorjl^HHl have been rigorously 
justified. ^SEZMl A rigorous proof of the leading order of 
the ground state energy per particle in the low density 
limit was given by Lieb and Yngvason.l^SESlxhe next or- 
der correction term was recently studied in a certain scal- 
ing limit.*" Bounds of the free energy at positive temper- 
ature were given iiP^ Cluster expansions give informa- 
tions on the phase without Bos e-Eins tein condensation, 
for repulsive or stable potentials.l^^IISl Recently there has 
been interest in Feynman cycles which should be related 
to Bose-Einstein condensation.*"' The conditions ([2]) and 
([3]) guarantee the absence of infinite cycles. This follows 
from the considerations here, and from the proof that all 
cycles are finite when the chemical potential is negative.*^ 

The following theorem gives bounds on the density 
p{z). Recall the function gr defined in ([6]). Let us define 
the following small parameter a(/3), which is associated 
with the scattering length a: for rf = 2, 

S(/3) - (1 Ho?IP)\ - 2 In I Ha/^) |) + | Ha'/PT^- 



and for d = 3, 



,1/2 



Theorem III.l. Let us assume that ln(a/v75)|~^ > 
_Ro when d ~ 2, or that a\f]5 > Rq when d = 3. Then we 
have, for < z < 1, 



p{z) ^ p^^Hz) j^^^ {h,{z)a{l3) + 2^^'a{P/2)) , 

(13) 



where 

hd{z) = (25 +5.(z)).9._i(z)+2i+\9.(z) + g.(z)2. 

(14) 

Notice that p(z) ^ p''^\z); this is an immediate 
consequence of (Hi). For d = 2 we believe that for z 
close to 1 the lower bound is optimal up to terms of 
higher order in a(/3), while for d = 3 the prefactor is 
not optimal. This is based on the (yet unproved) as- 
sumption that the leading order correction tothe pres- 
sure is equal to -87ra(/3)p(°)(-z)^ for z < l.EnEIISl. Us- 
ing (Is]) and ( [T2| , this suggests that p{z) « p'^^'>{z) — 
4a(/3j(47r/?)^~^/2(^)ffrf/2-i(2:)- If this indeed holds as 
a lower bound, one can replace the constant 5.09 in ^ 
by 3.52, yielding a bound in agreement with Huang's 
prediction.!^ 

From Theorem |III.l I we can easily deduce the bounds 
and ( 10 1, which we have used in the previous section. 

z/(l 



we see 



Since gQ[z) = z/(l — z) and gi{z) — — ln(l 
that the function h2 is bounded by 

h2{z) < C(l-z)-i|ln(l-z)| 

for some constant C < 00, which implies ([s]). To obtain 
(10 1, note that the function .93^(2) converges to C(3/2) 
1. Using the bound ([TT]) we see that hs is less 



as z 
than 



h,{z) ^ (23/2 -,C(i))y^+ 2^(1) + c(i)^ 



We are left with the proof of Theorem III.l In Sec 



tion |IV| we use the Feynman-Kac representation of the 
Bose gas to obtain bounds on the density. These bounds 
are expressed in terms of integrals of the difference be- 
tween the heat kernel of the Laplacian with and without 
potential. Section |V] deals with bounds of these inte- 
grals. It contains a novel variational principle for inte- 
grals over heat kernel differences (Lemma V.l[ ), which 
allows to bound these in terms of the scattering length of 
the interaction potential. Theorem |III. 1 1 then follows di- 
rectly from Proposition |IV.2| and from Lemmas |V.2| and 

EH 



IV. FEYNMAN-KAC REPRESENTATION OF 
THE INTERACTING BOSE GAS 

From now on we shall work in arbitrary dimension d ^ 
1. Let W* „ denote the Wiener measure for the Brownian 
bridge from x to y in time t; the normalization is chosen 
so that 



dWlyito) = (27ri)-'^/2 ^-\x^y\V2t ^ ^^(3. _ y). 



The integral kernel of e^'^^ — e'''^^"'^^ will be denoted 
by K{x, y). By the Feynman-Kac formula, it can be ex- 
pressed as 

K{x,y) = /(I - e--^^o'u{u>{s))As^^^Ap^^^^^ (^5) 
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Let us introduce the interaction U{ljj,lo') between two 
paths w and w': [0,2/3] W^. Namely, 



U{lo,uj') 



2/3 



U{uj{s)-uj'{s))ds. 



The following identity, which will prove useful in the se- 
quel, is obtained by changing to center-of-mass and rela- 
tive coordinates. 

Lemma IV. 1. For any x.y^x' ,y' e , 



2'^TT4p{x -~y + x' ~ y')K{x - x' ,y - y'). 



Proof. The difference w — a;' of two Brownian bridges is 
a Brownian bridge with double variance. Precisely, we 
have 



■^2f3{x - y) 



J Tr2/3{x' - y') 



4/3 



^0 



'Kij3{x- x' -y + y') 

By the parallelogram identity, 

'^2i3{x - y)T^2p{x' - y') 
TTip{x -x' -y + y') 



e 2 



-U(u,uj') ■ 



2'^nip{x~y + x' -y'). 



The result then follows from ( 15 I 



□ 



We also use the Feynman-Kac formula for the canoni- 
cal partition function. Namely, 

T.e-^H^-- = 4, E / dXi---dXN 



N 



exp 



1 ^ i<j < N 



(16) 



Here, Sn is the set of permutations of iV elements; Xj^{uj) 
is equal to one if u;(s) G A for all ^ s ^ 2(3, and 
it is zero otherwise. Eq. ( 16 1 makes sense for general 
measurable functions U : M U {oo}. In particular, 

we can consider the case of the hard-core potential of 
radius a. An introduction to the Feynman-Kac formula 
in the context of bosonic quantum systems can be found 
in Ginibre's survey."*^ 

We now rewrite the grand-canonical partition function 
in terms of winding loops. Let be the set of continuous 
paths [0, 2pk] M'' that are closed. Its elements are 
denoted by a; = {x, k, uj), with a; S M'^ the starting point 
and k the winding number; we have w(0) = Lo{2(3k) = x. 



For ^ £ ^ /c — l,we also let uj^ denote the £-th leg of 

ujt{s) = uo{2pi+ s), 
with ^ s ^ 213. We consider the measure /i given by 

k 

d/i(a;) = ^da;XAHdW^2'3fe(^)g-y(c.) 

Here, V{u}) is a self-interaction term that is defined below 
in Eq. ( [l7| . Let = Ufc ^ i^k', the measure jj, above nat- 
urally extends to a measure on fi. The grand-canonical 
partition function can then be written aa^ 

Z = — [ d^(a;i) • • •d^(a;„) 
xexp|- ^ 

1 z<j 5C n 

The self-interaction V{u)) and the 2-path interaction 
V{ijt},ijt3') are given by 



^ l<m ^ fc-1 



(17) 



We shall denote Vij = V{u}i,u}j) for short. Using (12 1 



one obtains an expression for the grand-canonical density, 
namely 



X j dM(w2) ■ ■ ■ j d/i(w„) e- ^'<^ . (18) 



From the representation ( 18 1 it is easy to see that 
p{z) ^ p'^'^\z). One uses the positivity of Vij to bound 
El ^<J ^ N E2 ^ ^<^, ^ N V'jj • For fixcd Wi, the 

integration over with j ^ 2 then yields Z, and hence 

p(z) ^ ^/ dM(a;i)fci s=: p(°)(2)- 

The last inequality follows since the self-interaction Viyi) 
is also positive. 

In the following proposition we shall derive a lower 
bound on p{z). We use the function defined in (14 1, 
as well as the integral kernel K{x,y) in (15 1. 



Proposition IV.2. For d ^ 1, (3 > and < z < 1, 

we have the lower bound 

p{z) ^ p(0)(z) - {ha{z) j K{x,y)dxdy 

+ j [K{x,x)+K{x,-x)]dx^ 

for any hounded (and measurable) A C M''. 
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Proof. Isolating the interactions between the first path 
and the others, we can bound exp{— X)i ^ i<j ^ n ^ij} 
from below as 



n 



J=2 



2 ^ *:</ n 



:5: 



l-X:(l-e-^«)]exp{- ^ Vm}. (19) 

i=2 2 ^ k<l ^ n 



We use this lower bound in Eq. ( 18 1. The first term (the 
1 in square brackets in (19)) is then |A|~^ J dfi{u>i)ki, 
since the integration over Ljj with j ^ 2 yields exactly 
Z. For the remaining terms (the sum over j), we also use 
the fact that the potential is repulsive so as to drop the 
interactions between ujj and the other loops in the last 
term in (19 1 for a lower bound. We conclude that 



Piz) > 1^ / d/i(a;)fc 



(20) 



1^ y'dM(a;i)fci I dM(a;2)(l- e-^-) 



In a similar fashion to ^T9[ . we have 



-y(a;) 



^1- (i-^ 

^ i<rn ^ k~l 

fcl-1 fc2-l 

g-V(wi,W2) ^ _ ^ _ g-C/(c^l.ti,"2,<!2) ^ ^ 

^1=0 £2=0 

Here, w^^f denotes the £-th leg of the path uji. We insert 
these inequalities into ( pO| , and obtain 

JO) 



-C/(cJ^,LJ„i) 



)> 



p(z) ^ p^^Hz)-A-i3, 



with 



^ £<m ^ A;-l 



-U{uJi,LJm) 



B^^J dK'^i)kl J d/i(a;2)fc2(l - e" 



■C'('^l,l,'^2,l) 



We also used xa ^ 1 to drop the restriction that paths 
stay inside A. Notice that only the first legs of ui and 
UJ2 interact in B; this is correct because we multiplied by 
kik2- 

We decompose the terms in A as ^1+^12+^3 according 
to the distance between interacting legs. Namely, the 
term fc = 2 in A is equal to 



z 



Using Lemma [IV. 1[ we get 



K{x, —x)dx. 



The terms with A: ^ 3 and two consecutive interacting 
legs are 

M = ^J^^dx.dx^dx:, I dW^^^,,^{u;^)dW^^^,^{uj2) 

Til,. 1 \ (fc + 2)z'''+2 |xi-^3l" 

^ (47r/3fc)'i/2 



Using Lemma |IV.1| and bounding the exponentials by 1 , 
we get 



A2 



2d/2^2 



gd/2-i{z) +2gd/2{z) K{x,y)dxdy. 



The terms where no consecutive legs interact are 
1 



A3 = —T / dxida:2dx3dx4 



E 



(fcl + fc2 + 2)^'^-i+'=2+2 



4,0k2 



, , ^ ^ (4^/3fci)'*/2(4^/3fc2)'^/2 

ki .k2 ^ i 

Then 

2''/2-l f |x-„-2.|2 

^3 I dxdydze s/j K{x,y) 



E 



(fci+fc2+2)z'^i+'^-2+2 



kiM > 1 



(fciA:2)'^/2 

-^gAz)\gi^{z)+gAz)\ I K{x,y)dxdy. 



(4^/3) 



We now decompose the terms in B as Bi + B2 + B3 
according to the winding numbers of toi and lu2. The 
term Bi involves two paths of winding numbers 1, and 
with the aid of Lemma IIV.II we find 



Next, B2 involves a path of winding number 1 and an- 
other path of higher winding number. Dropping the self- 
interaction terms yields the upper bound 



s$ 



2d/2 



X E l^d/2 



dxdy e K{x, y) 



{k + 2)z^+^ 



k ^ 1 
2 ^2 



2'^'-z 
(47r/3)' 



-i[9d/2-i{z) + 2gd/2{z)\ / K{x,y)dxdy 
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Finally, B3 involves paths with winding numbers higher 
than 2. We have 



2rf/2 
(4^ 



E 



_ |:c-y-2z|^ 

dxdydze K{x,y) 



(fcl + l)z'='-i+'=2+2 



/Ci,fc2 ^ 1 



(fclfc2)'^/2 

:^gdz)[9^_i{z) + gAz)] j K{x,y)dxdy. 



Collecting the bounds on Ai, A2, A3, Bi, B2, B3, we get 
the lower bound of Proposition |IV.2[ □ 



V. SCATTERING ESTIMATES 

As before, let U{x) ^ be radial and supported on 
the set {a; : \x\ ^ Rq}. Let a be the scattering length of 
U. We consider the Hilbert space L^(IR'') and the inte- 
gral kernel K{x, y) of the operator e^f^^ — ePi'^^-^) . 
It follows from the Feynman-Kac representation that 
K{x,y) ^ 0, see Eq. (151 in the previous section. 

We introduce 



1 

8^ 



K{x, y)dxdy . 



(21) 



We shall see below that, for d = 3, a{f3) is a good approx- 
imation to the scattering length. In fact, a ^ ^ 
flo, with ao the first order Born approximation to a. In 
two dimensions a{(3) is dimensionless and its relation to 
the scattering length is a{f3) « | ln(a^//3)|~^ for large f3. 
For i > 0, we also introduce the function 

1 „-t 

fit) = t 



t-1 



Lemma V.l. We have 



= — inf £^{^), 
Sir iieH^W) 



(22) 



where 



(2|V^/'(2;)p + U{x)\l ~ ijj{x)\^) dx 
^(^|/(/?(-2A + [/))|V). 



Note that / is monotone decreasing, with 1 ^ f{t) ^ 
2 for alH > 0. From monotonicity it follows immediately 
that a(/3) is monotone decreasing in f3. Moreover, for 
d = 3 it is not hard to see that lim^_»oo a(/3) — cl- For 
any d, lim^^o o,{P) = (Stt)"^ / U{x)dx. (This is also true 
when J U{x)dx = 00.) 

Proof. We first consider the case when U is bounded. 



With the aid of the Duhamel formula we have 



g2(/3-t)A^g2*A 



(3 i-t 



Hence 



32(/3-*)^C/e^(2A-C/){/e2(*--)Adsdi, 



a{P)^— j U{x){l-i^p{x))dx, (23) 



where ipf3{x) — {LpU){x), with 



Lp^ / (l-s//3)e^(2A-c/)d, 



The functional Sfj^ip) has a quadratic and a linear part 
in ^p, and it is not hard to see that the unique minimizer 
satisfies 



[-2A + U+y{P{-2A + U))]ij^U. (24) 



Since 



-"'ds 



it follows that tp — ipp, i.e. "0/3 = LnU is the unique 
minimizer of £p . After multiplying ( 24 1 by ipp and inte- 
grating we see that £i3{ipp) — J C/(a;)(T^-0^(a;))da; which, 
because of ( 23 1 , imphes (p2|) . 



Finally, the case of unbounded U can be dealt with 
using monotone convergence. If we replace U by Us{x) — 
min{J7(x), s} then the kernel K{x, y) corresponding to Us 
is monotone increasing in s. We can apply the argument 
above to and take the limit s — *■ 00 at the end. The 
convergence of a(/3) is guaranteed by monotonicity. □ 



The variational principle of Lemma |V.1| is convenient 
for obtaining an upper bound on a{[3). 

Lemma V.2. For d ^ 2 and ^/j3\ln{a/y/^)\-'^ > Rq, 



aid) < 



|ln(a2//3)|-21n|ln(a/V^)| " |ln(a2//3)|2 



For d ^ i and ay/J] > R\ 



(25) 



1 - (a/3i-'^/2) 



l-d/2Nl/(rf-l) 



(a/3 



1-^/2^,1/(^-1) 



(26) 



Note that the prefactor in ( |26| is equal to 1 for d — 3. 
Lemma |V.2| is the only place where the finiteness of the 
range Rq of U is being used. Appropriate upper bounds 
on a(/3) can also be obtained without this assumption. 
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and hence our main results generalize to repulsive inter- 
action potentials with infinite range (but finite scattering 
length). For simplicity, we shall not pursue this general- 
ization here. 



Proof. Let R > Rq, and let ■0oo be the minimizer of 



\x\ ^ R 



{2\V4>\^ + U\l-Tp\^)dx 



(27) 



subject to the bound ary condition tfjix) = for = R. 
It can be showip212U that there exists a unique minimizer 
for this problem, which satisfies ^ V'oo ^ 1 and 



in the region Rq ^ 
of (27 1 is given by 



ln(k|/a) 
ln{R/a) 
l-a\x\^-' 



for d = 2 
for d ^ 3 



^ R. Moreover, the minimum 



4tt 



p _ I ln(R/a) 



r(d/2)(l-a_R2 



for d = 2 
for d ^ 3. 



To obtain an upper bound on a(/3), we use the varia- 
tional principle (22 1 with ili{x) = ipooix) for |a;| ^ R, 



a.ndTp{x) — OfoT\x\ ^ R. Using |i/)oo| ^ land / ^ 2, 
we obtain the bound 



a(/3) sc: 



Er. ^ CFdR"^ 



where cr^ = 7r''/^/r(l + d/2) denotes the volume of the 
unit ball in M'*. The choice i? = 7/3 \ln{y/]3 / aVr^ for 
d = 2 and i? = (a^//3)l/('*-l) for d ^ 3 yields ^ and 

pel. □ 



For our lower bound on the density in Proposition |I V . 2"| 
we need a bound on two more integrals of the kernel K. 
Since they appear only in terms of higher order, a rough 
bound will do. 



Lemma V.3. Let 



a'(/3) = (87r/3)''/2-i / K{x,x)dx 



a" [13) = (87r/3)'^/2"i / K{x,~x)dx 



(28) 



The 



c{a'(/3),a"(/3)} < 2'''^a{l3/2) . 



For d = 3, it can be shown that both «'(/?) and a" {(3) 
converge to a as /3 ^ cxd, but we do not need this here. 

Proof. Using the semi-group property of the heat kernel 
we can write 



K{x,z)^ / (e'^^(x,2;)e'3^(2;,z) 

- e''(^-5^) (x,y)e''(^-5^) {y.z))dy. 

Since ab^cd ^ a(b— d) + b{a — c) ior a ^ c and 5 ^ d, 
K(x^ z) is bounded above by 



{x,y)(^e''^iy,z)- e^<^^-i-^^ {y,z))dy 
'^{y,z)[e^^{x,y)-ef'^^-'^^^ {x,y))dy. 



Using the bound e^'^ {x,y) < (47r/3) '^/^ the claim fol- 
lows easily. □ 



VI. CONCLUSION 

We have given rigorous upper bounds on the critical 
temperature for two- and three-dimensional Bose gases 
with repulsive two-body interactions. In two dimensions, 
our bound agrees to leading order in a^p with the ex- 
pected critical temperature for superfluidity. In three 
dimensions, our bound shows that the critical tempera- 
ture is not greater than the one for the ideal gas plus a 
constant times \fap^ . 

Our bounds are based on the observation that the one- 
particle reduced density matrix decays exponentially if 
the fugacity z satisfies z <\. What is needed are lower 
bounds on the particle density in the grand canonical en- 
semble. The Feynman-Kac path integral representation 
allows us to get bounds in terms of certain integral ker- 
nels which, in turn, can be estimated by the scattering 
length of the interaction potential using a suitable vari- 
ational principle. 

Acknowledgments: It is a pleasure to thank Rupert Prank 
and Elliott Licb for many stimulating discussions, and Markus 
Holzmann for helpful comments on the physics literature. D.U. 
is grateful for the hospitality of ETH Ziirich, the Center of The- 
oretical Studies of Prague, and the University of Arizona, where 
parts of this project were carried forward. Partial support by the 
US National Science Foundation grants PHY-0652356 (R.S) and 
DMS-0601075 (D.U.) is gratefully acknowledged. 



* Electronic address: rseirinK@princeton.edu ^ K. B. Davis, M. O. Mewes, M. R. Andrews, N.J. Van- 

t Electronic address: daniel@ueltschi.org druten, D. S. Durfee, D.M. Kurn, W. Ketterle, Phys. Rev. 

^ M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wie- Lett. 75, 3969-3973 (1995) 

man, E. A. Cornell, Science 269, 198-201(1995) ^ P. C. Hohenberg, Phys. Rev. 158, 383 (1967) 



10 



* V. N. Popov, Functional Integrals and Collective Excita- A. Siito, J. Phys. A 26, 4689-4710 (1993) 
tions (Cambridge University Press, 1987) "^^ D. Ueltschi, J. Math. Phys. 47, 123303 (2006) 

^ D. S. Fisher, P. C. Hohenberg, Phys. Rev. B 37, 4936-4943 K. Huang, Statistical Mechanics, 2nd ed. (Wiley, New 

(1988) York, 1987) 

" S. Pilati, S. Giorgini, N. Prokof'ev, Phys. Rev. Lett. 100, 

140405 (2008) 

R.P. Feynman, Phys. Rev. 91, 1291-1301 (1953) 

* T. D. Lee, C. N. Yang, Phys. Rev. 112, 1419-1429 (1958) 
^ A. E. Glassgold, A. N. Kaufman, K. M. Watson, Phys. 

Rev. 120, 660-674 (1960) 

K. Huang, in "Studies in Statistical Mechanics" , Vol. II, 
1-106 (North-Holland, 1964) 

A. Fetter, J.D. Walecka, Quantum Theory of Many- 
Particle Systems, Sec. 28 (McGraw-Hill, 1971) 

T. Toyoda, Ann. Phys. (N.Y.) 141, 154 (1982) 
" H. T. C. Stoof, Phys. Rev. A 45, 8398 (1992) 
^'^ M. Bijlsma, H.T. C. Stoof, Phys. Rev. A 54, 5085 (1996) 

P. Griiter, D. Ceperley, F. Laloe, Phys. Rev. Lett. 79, 3549 

(1997) 

" M. Holzmann, P. Griiter, F. Laloe, Eur. Phys. J. B 10, 
739-760 (1999) 

K. Huang, Phys. Rev. Lett. 83, 3770-3771 (1999) 

M. Holzmann, W. Krauth, Phys. Rev. Lett. 83, 2687-2690 

(1999) 

G. Baym, J. -P. Blaizot, M. Holzmann, F. Laloe, D. Vau- 
therin, Phys. Rev. Lett. 83, 1703-1706 (1999) 

^° J. Reppy, B. Crooker, B. Hebral, A. Corwin, J. He, G. 

Zassenhaus, Phys. Rev. Lett. 84, 2060-2063 (2000) 

V. A. Kashurnikov, N. V. Prokof'ev, B. V. Svistunov, Phys. 

Rev. Lett. 87, 120402 (2001) 
2^ P. Arnold, G. Moore, Phys. Rev. Lett. 87, 120401 (2001) 

H. Kleinert, Mod. Phys. Lett. B 17, 1011-1020 (2003) 

B. Kastening, Phys. Rev. A 69, 043613 (2004) 

K. Nho, D.P. Landau, Phys. Rev. A 70, 053614 (2004) 
'^^ G. Baym, J.-P. Blaizot, M. Holzmann, F. Laloe, D. Vau- 
therin, Eur. Phys. J. B 24, 107-124 (2001) 
J.-P. Blaizot, arxiv.org/0801.0009 

I. Bloch, J. DaUbard, W. Zwerger, Rev. Mod. Phys. 80, 
885-964 (2008) 

O. Bratteli, D. W. Robinson, Operator Algebras and Quan- 
tum Statistical Mechanics, Vol. 2, 2nd ed. (Springer, 1996) 

^° E.H. Lieb, J. Yngvason, J. Stat. Phys. 103, 509 (2001) 
E.H. Lieb, R. Seiringer, J. P. Solovej, J. Yngvason, The 
mathematics of the Bose gas and its condensation, Ober- 
wolfach Seminars (Birkhauser, 2005); cond-mat/0610117^ 

3^ F. J. Dyson, E. H. Lieb, B. Simon, J. Stat. Phys. 18, 335- 
383 (1978) 

T. Kennedy, E. H. Lieb, B. S. Shastry, Phys. Rev. Lett. 61, 
2582-2584 (1988) 

G. Roepstorff, J. Stat. Phys. 18, 191-206 (1978) 

V. A. Zagrebnov, J.-B. Bru, Phys. Rep. 350, 291-434 (2001) 
'■^^ J. Ginibre, Commun. Math. Phys. 8, 26-51 (1968) 

E.H. Lieb, R. Seiringer, J. Yngvason, Phys. Rev. Lett. 94, 

080401 (2005) 
^® A. Siito, Phys. Rev. Lett. 94, 080402 (2005) 

E.H. Lieb, J. Yngvason, Phys. Rev. Lett. 80, 2504-2507 

(1998) 

■'o A. Giuliani, R. Seiringer, J. Stat. Phys. 135, 915 (2009) 
■'^ R. Seiringer, Commun. Math. Phys. 279, 595-636 (2008) 
'^^ J. Ginibre, Some applications of functional integration in 

statistical mechanics, in "Mecanique statistique et theorie 

quantique des champs", Les Houches 1970, 327-427 (1971) 

S. Poghosyan, D. Ueltschi, J. Math. Phys. 50, 053509 

(2009) 



